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Abstract 

Finite size effects in the equilibrium phase space density distribution function 
are taken into account for alculations of the relaxation of collective motion 
in finite nuclei. Memory effects in the collision integral and the diffusivity 
and the quantum oscillations of the equilibrium distribution function in mo- 
mentum space are considered. It is shown that a smooth diffuse (Fermi-type) 
equilibrium distribution function leads to a spurious contribution to the relax- 
ation time. The residual quantum oscillations of the equilibrium distribution 
function eliminates the spurious contribution. It ensures the disappearance of 
the gain and loss terms in the collision integral in the ground state of the sys- 
tem and strongly reduces the internal collisional width of the isoscalar giant 
quadrupole resonances. 
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I. INTRODUCTION 



The relaxation of nuclear collective motion toward thermal equilibrium have been de- 
scribed in great detail within the framework of the kinetic theory, taking into account the 
collision integral In this theory, the damping of the collective motion appears due 

to the interparticle collisions on the dynamically deformed Fermi surface. It has already 
been shown by Landau [TO, IT] that both temperature and memory effects are extremely 



important for successful applications of kinetic theory to the relaxation processes in a Fermi 
system. However, only very little attention has been paid to the study of the peculiarities of 
the collision integral in a finite Fermi system caused by particle reflections on the boundary. 

In kinetic theory, the collision integral depends crucially on the phase-space distribution 
function f(f,p). The main aim of this paper is to apply the quantum Wigner phase-space 



distribution function, also known as the Wigner transform |]12| , to the evaluation of the 
collision integral in a finite Fermi system. The Wigner distribution function (WDF) is 
defined as the Fourier transform of the one-body density matrix over the relative coordinates. 



It possesses several nice properties pT3| , p!4| which justify its interpretation as the quantum 



mechanical analog of the classical phase-space distribution function. The WDF is useful in 



providing a reformulation of quantum machanics in terms of classical concepts fl!l^ , [T6|| and 
a good starting point for semiclassical approximations JT7 ] . 



Traditionally, the equilibrium phase space density distribution function f e q{r,p) in the 
collision integral is replaced by the one, feq^F^p), taken in Thomas- Fermi approximation 

feq,TF(r,P)=0(^-E(f,p)), (1) 

where 8(x) is the step function and E(f,p) is the classical single-particle energy. This is 
reasonable for an infinite Fermi system. In the case of a finite Fermi system, the quantum 
distribution function f e q(r,p) fluctuates strongly and contains the diffusivity of the Fermi 



surface even at zero temperature [|18|-f23|j . Both these features are due to particle reflections 
on potential walls. 
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The diffusivity of the Fermi surface in momentum space increases in the vicinity of the 
nuclear surface and thus enhances the effective particle scattering there because of the 
decrease of the Pauli blocking effect. It can be shown || that the use of a simple Fermi 
distribution function 

feqA^P) = (l + eX P 

with r-dependent diffusivity parameter a(r), instead of the ^-function of Eq. ([[]), for the 
equilibrium distribution function f eq {^,p) in the collision integral, leads to a significant 
enhancement of the damping of the nuclear giant multipole resonances. However, there is 
a conceptual disadvantage for the application of the Fermi distribution function f eq ,F(r,p) 
with a(r) 7^ to the collision integral. Namely, with this function the gain and loss terms 
in the collision integral are each nonzero for the ground state of the system, where the 
probability current should be absent by definition. We show in this work that in order to 
overcome this difficulty the smooth quantum distribution function f eq (r,p) should be used 
for feqi^iP) m the collision integral in the kinetic Landau- Vlasov equation. In contrast to the 
Fermi distribution function f eq ,F(r,p), the smooth quantum distribution function f eq (r,p) 
contains the residual oscillations [0,^] ensuring the above-mentioned condition for the 
disappearance of the gain and loss terms in the collision integral for the ground state and 
reducing the internal collisional width of the giant multipole resonances. 

In this paper we pay attention mainly to finite size and memory effects in the relaxation 
processes in finite Fermi systems. In Sec. II we obtain a general expression for the width 
of the giant multipole resonances at zero temperature of the nucleus starting from the 
collisional Landau- Vlasov equation. In Sec. Ill we study the influence of the memory effects 
and the diffusivity of the Fermi surface on the relaxation time. We use the smeared-out 
Wigner distribution function f eq {r,p) of the three-dimensional harmonic oscillator and of 
the Woods-Saxon potential for the calculations of the collision integral. A summary and 
conclusions are given in Sec. IV. 



E(f,p) — X 
a(r) 



\ ~ 1 



(2) 
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II. DAMPING OF COLLECTIVE EXCITATIONS IN KINETIC THEORY 



The kinetic equation for a small variation 5f(f,p,t) of the distribution function can be 
transformed to a set (infinite) of equations for the moments of Sf(f,p, t) in p space, namely, 
the local particle density 5p, velocity field u, pressure tensor 7r a p, etc., see O,^!]. The 
first-order moment of the kinetic equation has the form of the Euler-Navier-Stokes equation 
and is given by @.^| 

d d (5 2 £\ d 

mp eq -u a + Peq — jjp+ —n ua = 0, . (3) 

Here and in the following expressions, repeated greek indices a, (3, v = 1,2,3 are to be 
understood as summed over. The variation of the local particle density p and the velocity 
field u a in Eq. are defined by 

Sp = 9 IjS^ Sf > * = 9 Jj^wi 6f - (4) 

The quantity 7r a/ g is the deviation of the pressure tensor from its equilibrium part P eq due 
to the Fermi-surface deformation 

71-0/3 = ml (2i!h) 3 ^ Pa ~ mUa ^ Pp ~ mu p) 5 f ( 5 ) 
and the equilibrium pressure, P eq , is given by 

_ r dp p 2 ^2 
^- g J {2Tihf2m hq - 3 kmi { ' 

where S^m is the kinetic energy density. The internal energy density £ in Eq. ([3D contains 
both kinetic and potential energy densities: £ = £y. in + £ pot , where £ pot is the potential 
energy density. 

Eq. ([|) is not closed because it contains the pressure tensor 7r a/3 given by the second- 
order moment of the distribution function 5f(r,p,t), Eq. ([5]). We will follow the nuclear 
fluid dynamic approach of Refs. 0,[||6| and take into account the dynamic Fermi-surface 
distortions up to multipolarity 1 = 2. The second p moment of the kinetic equation leads then 
to a closed differential equation for the pressure tensor 7r Q/ g. Namely (see Refs. []3|,|6|,|8|J^,p6[]), 



9 „ fdu a f dup 2 du a \ ir a/3 



The local relaxation time r 2 in Eq. (0) is caused by the interparticle scattering on the 
deformed Fermi surface: 

dpp 2 Y 20 5 St 



1 



/ A mr? 



dpp Y 20 Sf 



(8) 



Here 5 St = SSt(f,p,t) is a collision integral linearized in 5f. In the case of small eigenvi- 
brations with eigenfrequency uj = ujq + iT/2h, where u>o and T are real, the collision integral 
can be transformed, taking into account also memory effects, as []5|j27f] 



8St(f,p, t) 



g dp 2 dp 3 dp 



'4 
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j=l °Jj eQ 1 

where W({pj}) is the probability of the scattering of nucleons near the Fermi surface, g = 4 
is the spin-isospin degeneracy factor and 

Q = (l - /0(i - f 2 )hh - /i/ 2 (i - / 3 )(i - / 4 ); 

AE = E ± + E 2 — E 3 — E^ Ap = p 1 +p 2 -p 3 - p 4 ; 

Ej = p^/2m + V(rj) is the classical single-particle energy, pi = p and V(r) is the nuclear 
mean field. 

We will follow the arguments of the Fermi-liquid theory of Ref. pi and assume that the 



dynamical component of the distribution function Sf(f,p,t) has the form 

6f(r,p,t) = -^v(r,p,t), (10) 

where u(r,p,t) are unknown functions. The functions u(f,p,t) depend on the orientation p 
only because of the sharp energy dependence of the factor df eq /dE in Eq. flID"D, which is 
localized at the Fermi momentum Pf{t)- We point out that the smooth region of the equi- 
librium distribution function f e q{r,p) in momentum space appears in a classical forbidden 
region at E < V(r). However, this region is absent in the space integral in Eq. Thus, 



v(r,p,t) « u(r,p F (r),p,t) =J2 u i™(r,t)Y lm (p). (11) 



Im 



An exact evaluation of the nine- dimension integral in p space in Eq. @ is a very 
complicated problem. We will follow the Abrikosov-Khalatnikov method |24| improved in 
Ref. |28fl . Let us assume that the scattering probability W({f>j}) in Eq. (|§D depends on two 



scattering angles 9 and only, (see Ref. p4|), where # is the angle between pi and p 2 , and 
is the angle between the planes formed by (pi,f>2) and (^3,^4), i.e., 

W(iPj}) « = Mr), 0, 0})- (12) 

To evaluate the collision integral Eq. @, we will use the Abrikosov-Khalatnikov transfor- 



mation 24,28 



r dp 2 dp 3 dp 4 m 3 r dfid0 2 

y w ( pl ( -° 272W y ^m dE2 ( - } ' (i3) 

where <if2 = sin 8d8d<j) and 02 is the azimuthal angle of the momentum p 2 in the coordinate 
system with the z axes along p\. We point out that the angle varies only from to 7r 
because the particles are indistinguishable. 

Using Eq. (|H|), the collision integral Eq. (||) can be written in the form 



8St(r,p, t) 



3 

m 



16vr 4 ^ 



X 5>m(r,*)*bn(p) EWM)fl(cosej)) (/| +) + jj^) . (14) 
l,m j=l 

The symbol (...) denotes the averaging over angles of the relative momentum of the colliding 
particles 

(w(e,m(cose 3 )) = [dB-^j- r£w(e,m{<x*0j), (is) 

JO COS u / Z JO Z7T 

where cos9j = (pjPi), i.e., # 2 = and 

cos 3 = (cos(#/2)) 2 + (sin(#/2)) 2 cos 0, cos # 4 = (cos(fl/2)) 2 - (sin(#/2)) 2 cos 0, 
P;(cos#) is a Legendre polynomial and 

3 Jv eq OEj Sfj ^ 
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We now return to the dynamical equation fl3|). Let us introduce the displacement field 
x(r , t) related to the velocity field u(f, t) as 



->r-> \ 9 _._ . 

We look for the displacement field x(f, t) in the following separable form: 



X(r,t) =P{t)v(r), 



(16) 



(17) 



where (3{t) = j3 e luJt . Substituting Eqs. (|IIJ) and ([T7|) into Eq. @, one obtains the equation 
of motion for the collective variable @(t): 



1 

mp eq vj + -— (Imn^J $ + — (RelT^) (3 



-Peg 



_8_ 



'5 2 £\ d 



5p 2 I dr 

■ ' /eg 



13 = 0, 



(18) 



where we have used the following form for the traceless part 7r Q/ g(r, t) of the momentum flux 
tensor 



(19) 



Multiplying Eq. (0) by v a , summing over a and integrating over r space, we obtain the 
dispersion equation for the eigenfrequency uj\ 



- Buo 2 + iuA(u) + C{u) = 0. 



(20) 



Here, B is the hydrodynamical mass coefficient [2S] with respect to the collective variable 



B = m J drp eq v 2 . 



(21) 



The dissipative term A(u) and the stiffness coefficient C{uo) 





C + C'(lo) are given by 



• 4(w) = ^/ <i "'«*: (Imn; 



(22) 



and 
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c 



df 



> 2 , 



_d_ 



(peqVuJ 



C'(u) = I drv a — (RelT 
J or u 



(23) 
(24) 



We point out that the definition of the stiffness coefficient C coincides with the one for 
the stiffness coefficient in the traditional liquid drop model (LDM) for the nucleus. In 
contrast, the additional contribution C'(uj) to the stiffness coefficient is absent in the LDM 
and represents the influence of the dynamical Fermi-surface distortion on the conservative 
forces in the Fermi system. Finally, the dissipative term A(u) appears due to the interparticle 
scattering on the distorted Fermi surface. 

To evaluate the coefficients A{uS) and C'(u), we will use the third equation of motion 
(0). Let us rewrite Eq. (]?]) in the form 



where 



d 1 

dt H t 2 



du a dup 2 du a 

A-a/3 = H -On 



(25) 



drp ~^ dr. 



Ja/3 



drr- 



(26) 



Taking Eqs. fll6|) , C p~9| ) and (]25[) , one obtains 



/3ReIT^ + -/3Imn^ + -/3Ren^ + — /3Imn^ = -P eq $A a/3 . 



(27) 



Here. 



A, 



aj3 



dv n dv 



P 



--6. 



dVr 



drp ' dr a 3 a/3 'dr 



(28) 



Eq. (|27D can be represented as a set of equations of motion for both the real and the 
imaginary parts of oj: 



2h 



- TReU^ - 2HM^ + —ReU%g - TP eq k aP = 0, 

2h - 

2hw ReW, - rimU%g + —IraW, + 2huj P eq K ap = 0, 

r 2 



(29) 
(30) 
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where Im(u;) = T/(2h). 



In the case of small damped collective motion, we find from Eq. (p7[) 

Ren^ = -P eg i ^g- A^, (31) 



ImlT^ = -P et . ° 2 , 2 V (32) 



1 + (oo r 2 ) 2 
1 + (u; t 2 )' 

Finally, we have from Eqs. fl22|), (|§, (0) and (see also Ref. f2|), 

^-J^TTT^^' (33) 

In the same case of small damped collective motion, the dispersion equation (^) is 
transformed as 

r 

- B(ujI + i-uo) + iuj A(uj ) + C(u ) = 0. (35) 



Thus, 



= — g— , r = g— • ( 36 ) 



Using Eqs. (0), (0) and (|34|) we obtain the width T as 

uj A(uj ) 



hu - 



CM 



dfP eq A ai/ (dv a /dr u )uj T 2 /[l + (c^o^) 2 ] 

huj J — f . (37) 

C + / drP eq Ku(dv a /dr v )(u r 2 ) 2 /[l + (uj r 2 ) 2 } 



In the case of an isoscalar giant quadrupole resonance (GQR) one can assume P,^JT7| 

that the displacement field v{f) is the irrotational one and is given by v = (—x,—y,2z). 
Qu 

We have then A au —^- = 8. Furthermore, the LDM stiffness coefficient C gives a negligible 

or u 



contribution to the total value C(u) f3|,|171] and Eq. (|37| ) is transformed as 



(^oT2) 21 



/ dfP eq u T 2 /[l 

T = ^o-f • (38) 

dfP eq (u T 2 ) 2 /[l + (LU T 2 ) 2 } 



Finally, in the rare collision regime (ujqt ^> 1) Eq. (|3^) is reduced as 

P. 



r, - // / df'-^-j J drP eq . (39) 



III. RESULTS OF NUMERICAL CALCULATIONS 

To study the space distribution of the collective damping, we will introduce the local 
damping parameter £(r, uj ) related to the width T: 

T = J df%(r,u ), (40) 



where (see Eq. (|38p) 



£(r, u ) = ^o^^or 2 /[l + [u a r 2 f] ^ 



drP eq (co r 2 y/[l + (lu t 2 ) 2 ] 

is always a positive quantity. The local relaxation time t 2 = r 2 (r,uJo) can be obtained by 
substituting Eq. ([14]) into Eq. (§) and is given by 



1 m 3 4- 

" dW(9A)) ,oo ^ = 7 (r,wo), (42) 



r iE(E -v„r-% 



where 



roo 

R {±) = / dE 1 dE 2 dE 3 dE 4 (E 1 - V eq ) 3/2 S(AE ± hu ) 
Jv eq 

„(df eq ,i SQ df eqa 5Q ,N<9/eg, 3 5Q \ 

x k^ +C! «^ +(1+ci ^Wfs)' (43) 

and the coefficients c 2 and d 2 are given by 

c 2 = (WP 2 (cos^))/(W r ), d 2 = (3 W sin 4 - sin 2 0)/(PF). (44) 
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To evaluate the relaxation time fl4"2"|), we will study, first of all, the collision integral 
St eq (r,p) in the ground state of the system which is given by (see also Eq. (||)) 

St eq {m = J 9d ff^f 4 Q \j(AE) 6(Ap). (45) 
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We will introduce also the total gain and loss fluxes of the probability in the ground state. 
They are given, respectively, by 

Je q ^n = J dn ^^^f^ Wm}) [l-f eqA ][l-fe q Me q ,3fe qA 6(AE)5(Ap), (46) 



and 



Jewess = J Pl T\y ^ W({Pj}) [1 " /Jll " /eftj W W (47) 

One can see from both definitions, Eqs. (filf ) and (|47|) , that J eq ^ a in = Je q ,ioss and St eq (r,p) = 
0, as it should be for the equilibrium state of system. Moreover, in the case of the ground 
state of the system, both fluxes J eq , g ain and J eq ,ioss have to disappear separately. This is 
not the case, however, if the Fermi distribution function f eqF (see Eq. (|2|)) is used for the 
equilibrium distribution function f eq of the ground state of the finite Fermi system in Eqs. 
(f|6|) and (f47|). To avoid this disadvantage we will use the smeared-out quantum distribution 



function f eq in both Eqs. (^) and (0). To control the disappearance of the gain and loss 
fluxes, we will introduce the relative contribution q of the probability fluxes J eq , g ain or J eq j OS s 
evaluated with f eq = f eq to the corresponding values evaluated with f eq = f eq ,F, i-e., 

<? = Je q ,gain(\feq}) / Jeq,gain({fe q ,F}) = Je q ,loss ({fe q })/J eq ,loU{fe q ,F})- (48) 

Below we will apply our approach to the isoscalar GQR. We will assume that the scat- 
tering probability W in Eqs. (^2|) and (f44|) is angle independent, i.e., di = 4/5 and C2 = 1/5, 
and the magnitude of W can be obtained from the nuclear matter estimate of the parameter 
a = 15 7r 2 h 5 /m 3 (W) = 9.2 MeV from Ref. ||. We will also use the following expression 
for the energy Tiuq of the isoscalar GQR, Tlojq = 6(M _1//3 MeV. The numerical calculations 
will be performed for both the spherical harmonic oscillator (HO) potential well and the 
Woods-Saxon (WS) potential. 

A. Spherical harmonic oscillator potential 

We will use the harmonic oscillator potential in the form 
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V eq (r) = ^mOV, (49) 



where Ml ~ 41A" 1 / 3 MeV. For " magic" nuclei in the absence of the spin-orbit interaction, 
the smeared-out quantum distribution function f e q{r,p) in a HO potential is given by Jl9] 



f eq (f,p) = f eq (e) = 8e- e £(-l) fc ^(2e)n fc . (50) 

k=0 

Here, L%(e) is the associated Laguerre polynomial and e = p 2 /mhQ + mQr 2 /h = 2E/hQ is 
the dimensionless energy parameter. The smooth occupation numbers n k are introduced as 

dx((x) + 

where ((x) is the averaging function chosen as 



n k = dx((x) + 2^ a2 M C2 M -i(Afc), (51) 

J— CO i 



C(x) = \e-*\ 



The second term in Eq. (|5i~D contains the so-called Strutinsky curvature corrections and is 
given by 

Ux) = K -^Le~ x H n (x), 

\/7T 



where H n (x) are the Hermite polynomials, a 2n = (— l) n /(2 2n n\), and = (Ep — E k )/^. 
The quantities Ep and are the Fermi energy and single-particle energies in the mean 
field V eq (r), respectively, and 7 is an averaging parameter. 



The results of numerical calculations (see also Ref. [I9[) of the smooth distribution func- 
tion f eq (e) of Eq. (|50D are shown in Fig. 1. The solid line 1 gives the behaviour of the smooth 
distribution function for the value of the smearing parameter 7 = 2.5 • hQ and 2M = 6, the 
order of the curvature correction polynomial in Eq. (0). We point out that these values 
of the smearing parameter and correction polynomial are localized in the so-called plateau 
region for the shell correction 5U to the binding energy, i.e., where 5U does not depend on 
7, Ref. [0. 

The smooth distribution function f eq (e) exhibits oscillations caused by particle reflections 
on the potential surface. The mean behaviour (i.e., without the oscillations) of f e q{t) can 
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be approximated by the Fermi function / e?) i?(r,p) of Eq. (0), where the chemical potential 
A is fixed by the condition of conservation of particle number A: 

dfdp 



A = 9 J 



:feqAr,P)- 



(2vrft) 3 ' 

The solid line 2 in Fig. 1 shows the behaviour of the Fermi distribution function f e q,F{r,p) 
of Eq. (0) with the parameter a taken from Eq. ([55]). The dashed line in Fig. 1 gives the 
simple Thomas- Fermi approximation Eq. ([I]). 

In Fig. 2 we show the smooth distribution functions f eq as functions of the dimensionless 
kinetic energy €ki n = p 2 /mTiVt for different distances r. We can see from this figure that the 
diffusivity parameter a is almost independent of the distance r. This fact is a feature of 
the HO potential well. Note that the analogous diffusivity parameter for the Woods-Saxon 
potential is a strongly r-dependent function increasing near the potential wall, see below 
and Ref. |2"2"|,|2"3"|| . Following Ref. P3f, the diffusivity parameter a of the quantum distribution 



function f eq in momentum space can be estimated using the expansion of f eq in Hermite 



polynomials. The result reads 





G 2 + 


\G 3 - 


2/3 


CL — \ 










. 2 . 





(52) 



The parameters G 2 and G3 depend on the mean field V(r) and are given by (in the lowest 



order of ftf 



G. = -*- 
4m 



zv: q (r) + K» 



-.2 r 



3 4m 



(V' (r)) 2 + 2- 
y eqK n 3m 



(53) 



where prime means an r derivative. For the HO potential one has G 2 = — 3(hQ) 2 /A and 
G3 ~ — (hQ) 3 X/2. The diffusivity parameter a can be also estimated from a fit of the smooth 
distribution function f eq and its derivative df eq / dE to the corresponding values of the Fermi 
distribution function Eq. (0) within some smeared out interval 7 < 7 near the Fermi energy 
Ep. Namely, one has the following estimate: 



7 Je f - 7 /2 df eq /dE 



(54) 
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Note that Eq. ( (54"|) gives an exact result for the diffusivity parameter a if the distribution 
function f eq coincides with the Fermi distribution function of Eq. (0). 

The solid lines 1 and 2 in Fig. 3 show the numerical results for the parameter a as a 
function of mass number A obtained using both expressions (|52|) and (0), respectively. A 
discrepancy between the results presented by curves 1 and 2 appears because the expression 
(|52|) takes into account the lowest orders of expansion of the distribution function f eq over 
the Hermite polynomials; see Ref. [ 23[| . We have also established the following simple A 



dependence for the diffusivity parameter a for the spherical HO potential well, obtained 
from Eq. @ with 7 = 7/2 = 1.25 HQ, 

a « 10.2A~ 1/6 MeV. (55) 



The result of the numerical calculation of the diffusivity parameter a Eq. fl55|) is shown in 
Fig. 3 by the dashed line. 

We will return now to the problem of evaluating the relaxation time, Eq. In Fig. 4 we 
have plotted the ratio q as obtained from Eq. ( f4"8| ) as a function of the mass number A (curve 
HO for the harmonic oscillator potential ([Sf)). We have used here the Fermi distribution 
function Eq. (Q) with parameter a from Eq. (^). As can be seen from Fig. 4, the effect of the 
quantum oscillations for the smeared-out distribution function f eq = f eq leads to an essential 
compensation of the gain and loss probability fluxes. A small nonzero contribution (about 10 
- 20 % in Fig. 4) remains in part because of the Abrikosov-Khalatnikov transformation Eq. 
(p~3|) used earlier, which implies a localization of the momentum pj near the Fermi surface 
in Eq. (|9|). In order to check this statement, we will use more general transformation of 
the momentum integrals in the collision integral, Eq. (|9|), with an arbitrary value of the 
momentum pj, suggested in Ref. 0: 

dp 2 dp 3 dp A -^r/ AP w n 

P1+P2 + 2pip 2 n - Jpj + p\ - 2p 1 p 2 n (...), (56) 



3 

m 



167r 4 r Pl 



— i — / dEodE-] 



where n = min (l, \Jplipi + p\ — pD/ViVt) ■ The solid line 1 in Fig. 5 gives the behaviour 
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of the ratio q as a function of the smoothing parameter 7 as obtained from Eq. (fi8| ) by 
applying the transformation Eq. ([56]) to the calculations of the probability fluxes, Eqs. 



( |i6D and (f4?D. A strong deviation of this result for q from the analogous one obtained 
using the Abrikosov-Khalatnikov procedure, Eq. ([OJ), (the solid line 2) appears for small 
magnitudes of smoothing parameter 7. This is because the smeared-out distribution function 
f eq oscillations grow with the decrease of the smoothing parameter 7 and the accuracy of 
the Abrikosov-Khalatnikov transformation Eq. (|l3| ) also decreases. 



The results of the numerical calculations of the local relaxation time (more precisely, the 



inverse value h/T2(r 1 uo)) of Eq. (f42|) for the isoscalar GQR in the nucleus with A = 224 
are shown in Fig. 6. We point out that the use of the Fermi distribution function Eq. (0) 
instead of f eq in Eq. (|42] ) (the solid curve 2 in Fig. 6) leads to much stronger damping than 
the analogous calculations with f eq from Eq. (0) (solid curve 1). 

To give a simple phenomenological prescription for the removal of the nonphysical prob- 
ability fluxes (|46| ) and ( f4"T|) in the ground state of a system, we will introduce the modified 
distribution function 

Feq = feq,F + V^feq- (57) 

Here, f eq ^ is the Fermi distribution function Eq. (0) with diffusivity parameter a from Eq. 
(H), and Af eq is given by 

A feq feq feq,Fi 

where f eq is the smooth distribution function of Eq. (|50D . A numerical calculation of the 
probability fluxes J e q,gain, Eq. (gfj), and J eq ,ioss, Eq. ((47]), with f eq replaced by T eq shows 
that both probability fluxes disappear at rj = 0.86 for the nucleus with A = 224. The 
corresponding distribution function T eq is shown in Fig. 1 as the dotted line. The local 
relaxation time from Eq. ( |42"D obtained with modified distribution function from Eq. ( |57|) 
is shown in Fig. 6 as solid curve 3. It should be remembered that the two last results 
(curves 1 and 3) are more correct in the sense of the compensation of the probability fluxes 
in the ground state of the system. Both these results show strong oscillations of the local 
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relaxation time within the nuclear volume and amplification of the damping in the surface 
region of the nucleus. This feature of the finite Fermi system arises due to the fact that 
the equilibrium distribution function f eq , Eq. (|50| ) fluctuates strongly and contains the 
diffusivity of the Fermi surface. It is interesting to note that, excluding the oscillations in 
the nuclear volume, the result noted as curve 3 in Fig. 6 with full compensation of the 
probability fluxes in the ground state is in good agreement with the one (^/t 2 )tf obtained 
using the simple Thomas-Fermi distribution function of Eq. ([!]) (the dashed line in Fig. 6 
with {h/T 2 ) T F = 3(^ ) 2 /(4vr 2 a)). 

We point out that the quantum calculations represented in Fig. 6 by curves 1 and 3 
have been done with a quite large smoothing parameter 7 = 2.5 • hQ. A decrease of 7 leads 
to negative values of the local relaxation time r 2 (r, u ) in some regions of r. The behaviour 
of the damping factor £(r, uq) is shown in Fig. 7. It is necessary to stress that nonzero 
damping in a cold Fermi system appears only because of the memory effects in the collision 
integral (|9]). 

The results of numerical calculations of the width T of the isoscalar GQR as a function 
of the mass number A are shown in Fig. 8. As can be seen from this figure, the smooth 
distribution function f eq from Eq. (^) leads to the contribution of the collisional relaxation 
to the isoscalar GQR width (solid line 1) which does not exceed 30-50% of the experimental 
values. This result is similar to the one obtained with the sharp Thomas-Fermi distribution 
function Eq. ([!]) (the dashed line) and agrees with the earlier calculation of the internal 
collisional width of the isoscalar GQR ||. In contrast, the analogous calculation with the 
smooth Fermi distribution function (Q) (solid line 2) largely overestimates the contribution 
of collisional damping. Note that the small nonzero probability fluxes appearing in J eq ,gain 
(f46|) and J eq ,ioss (0) evaluated with f eq from Eq. (|50| ) lead to a very small contribution to 
the final result for the width. This can be seen in Fig. 8 from a comparison of both curves 
1 and 3. 

B. Woods-Saxon potential 
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We have applied above the HO Wigner distribution function to the calculations of the 
collision integral in a finite Fermi system. This distribution function contains both important 
ingredients of influence of the multiple particle reflections from the potential surface: the 
diffusivity and the oscillations of the distribution function in momentum space. However, 
the realistic nuclear potential well has a finite depth, providing stronger surface effects than 
the ones in the HO mean field. We will give below an analysis for the case of the WS 
potential in the form 

Vws(r) = V /(l + exp[(r - R )/d}). (58) 

We have adopted the following parameters Vq = —44 MeV Rq = 1.27 A 1 / 3 and d = 0.67 fm. 

An exact quantum calculation of the equilibrium distribution function f eq (r,p) for the 
WS potential is a rather complicated problem. We will use the result of Ref. |Qj for the 
semiclassical expansion of f e q{r,f)) in the Hermite polynomials. This gives 



N 



n=0 ™ ! V O- 



fe q (r,P) = E ( -T 1 ) X=E - (59) 



Here <&(y) = [l+erf(y/ v2)]/2 is the normal distribution function, erf (y) is the error function, 
$(o) = <^ anc [ j nc l ex n a t $( n )((x — fJ,)/cr) denotes the x derivative of n order. We will 
take into account terms with n < 3 in the expansion fl59|) . The corresponding coefficients b n 
are given by 

b = 1; h = E(f,p) - n; b 2 = G 2 + bl- a 2 ; b 3 = G 3 + b\ + Zb^Gz - a 2 ), 

where G 2 and G3 are obtained from ((53|). The coefficients \x and a can be found by solving 
the following nonlinear equations: 

6 4 = b\ + 6bl(G 2 - a 2 ) + 46x^3 + 3a 2 (a 2 - 2G 2 ) = 0, 

65 = 15a\ - 10a 2 {bf + 3b x G 2 + G 3 ) -b\- l$b\G 2 - 106 2 G 3 = 0. 



Following Eq. (g), the mean behaviour of f eq (r,p), Eq. fl5T3|), can be approximated by 
the following Fermi function: 
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feq,ws(r,p) = 1 +exp 



E(r,p) — A 
a(r) 



(60) 



Here E(f,p) is the classical energy of the particle in the WS potential 



P 



E(r,p) = ^- + V ws (r) 
2m 



The local diffusivity a(r) is evaluated from Eq. 
approximation (V^(r) « V^{Rq) = 0) @ 



with taken in the so-called surface 



g 2 = -^i 

4m 



— V (r) 



4m 



where pp{r) = h(3n 2 p eq (r) /2) 1 / 3 is the Fermi momentum taken in the local density approx- 
imation. The results of numerical calculations of the parameter a(r) are plotted in Fig. 9 
for both HO and WS distribution functions. We can see that the parameter a(r) for the 
WS potential well is peaked in the vicinity of the nuclear surface. The distribution function 
f e q(r,p) smeared out over angles in p space is plotted in Fig. 10. It contains both diffusivity 
and oscillations in p space which depend on the distance r. In Fig. 4 we show the ratio q as 
defined by Eq. (H3) (curve WS). We have used there the Fermi distribution function, Eq. 



fl60D, instead of the smooth distribution function f eq ^ in Eq. (figf). One can see that the 
occurrence of the quantum oscillations in the smeared-out distribution function f eq leads to 
an essential compensation of the gain and loss probability fluxes. 

The local relaxation parameter (|42D and the damping factor (41) are shown in Figs. 11 
and 12, respectively. It can be seen that the collisional damping is more pronounced in the 
case of the WS potential than in the HO one shown in Figs. 6 and 7. We point out also that 
the occurrence of the quantum oscillations in the equilibrium distribution function strongly 
reduces the relaxation processes in the finite Fermi system. The total collisional width T of 
the isoscalar GQR as a function of the mass number A evaluated for WS potential is shown 
in Fig. 13. The final result (solid line 1) agrees with above calculation of the collisional 
width in the HO potential; see the dashed line in Fig. 8. 
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IV. SUMMARY AND CONCLUSIONS 



Starting from the collisional kinetic equation we have derived the isoscalar GQR width 
for a finite nucleus taking into account the memory effects and the peculiarities of the equilib- 
rium distribution function f eq caused by the multiple reflections of particles on the potential 
wall. The equilibrium distribution function contains both smooth and the oscillating com- 
ponents. The smooth component of f eq can be approximated by the Fermi distribution 
(0) with an r-dependent diffusivity parameter a(r) in momentum space. The diffusivity 
parameter a(r) is almost independent of the distance r in the case of the harmonic oscil- 
lator potential well and is a strongly r-dependent function in the case of the Woods-Saxon 
potential, increasing near the potential wall. We have demonstrated numerically that the 
smooth part of the equilibrium distribution function can be satisfactorily described by using 
the semiclassical expansion of f e q{r,p) on the Hermite polynomials; see Ref. 



It was shown that the general condition for the disappearance of the gain and loss 
probability fluxes in the ground state of system can be reached due to the occurrence of 
quantum oscillations in the distribution function in momentum space in the nuclear volume. 
The diffuse tail of the distribution function in momentum space leads to an increase of the 
collisional damping of the collective motion in the surface region of the nucleus and thus 
to an increase of the isoscalar GQR width. However this increase is strongly reduced due 
to the above-mentioned oscillations of the equilibrium distribution functions appearing in 
the collision integral Eq. (g). As a result the collisional width of the isoscalar GQR does 
not exceed 30-50% of the experimental value and agrees with the estimates of the width 
where the sharp Thomas- Fermi distribution function (i.e., in the absence of the diffusivity 
and quantum oscillations) is used. The collisional damping in a cold Fermi system arises 
only because of memory effects in the collision integral (|9]). 

To describe the experimental values of the multipole giant resonances additional con- 
tributions from other spreading sources, such as the fragmentation width in random phase 
approximation calculations or its representation through the one-body dissipation (see Ref. 
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H) have to be taken into account. We pointed out also that surface effects in the collisional 
damping are manifested more distinctly in the Woods-Saxon potential where the diffusivity 
parameter a(r) of the distribution function in momentum space is r dependent and increases 
within the surface region of the nucleus. 
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FIGURES 




FIG. 1. The equilibrium distribution function as a function of the dimensionless parameter e, 
for a nucleus with mass number A= 224, calculated for a spherical HO potential. Solid lines 1 
and 2 show the smooth distribution function obtained using an averaging procedure of Eq. (50) 
with 7 = 2.5 • HQ and the Fermi distribution function of Eq. (2) with parameter a from Eq. (55), 
respectively. The dashed line shows the Thomas-Fermi distribution function (1). The dotted line 
shows the distribution function of Eq. (57) with rj = 0.86 providing the disappearance of the 
probability fluxes in the ground state of the nucleus. 
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FIG. 2. The smooth distribution function of Eq. (50) as a function of the dimensionless kinetic 
energy eki n = p 2 /mfiQ, for a nucleus with mass number A= 224, calculated for a spherical HO 
potential with 7 = 2.5 • hfl. The different curves correspond to the distances r = 0, 3, and 6/m to 
the center of the nucleus. 
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FIG. 3. The diffusivity parameter of the smooth distribution function of Eq. (50) versus the 
mass number A, calculated for a spherical HO potential with 7 = 2.5 • HQ: For curve 1 we use Eq. 
(52), for curve 2 we use Eq. (54) and the dashed curve is obtained from the fitting formula of Eq. 
(55). 
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FIG. 4. The relative probability fluxes as given by Eq. (48) versus the mass number A for both 
HO and WS potentials. 
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FIG. 5. Dependence of the relative probability fluxes q on the averaging parameter 7 in units 
of hfl: For solid line 2 we use the Abrikosov-Khalatnikov transformation Eq. (13), and for solid 
line 1 we use the transformation Eq. (56), for the momentum integrals in the collision integral Eq. 
(9). The calculations were performed for the nucleus with A=224 in the HO potential well. 



28 




FIG. 6. The local damping parameter h/T2(r,u>o) for the nucleus with A = 224, calculated for 
a spherical HO potential. The different curves correspond to the different equilibrium distribution 
function in Eq. (42): For curve 1 we use the smooth distribution function of Eq. (50) with 
averaging parameter 7 = 2.5 • hQ, for curve 2 we use the Fermi distribution function (2) with a 
from Eq. (55), and for curve 3 we use the distribution function from Eq. (57). The dashed line 
shows the result obtained using the Thomas-Fermi distribution function of Eq. (1). 
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FIG. 7. The damping factor £(r, uo) as a function of the distance r. The notations is the same 
as in Fig. 6. 



30 



8 



> 



6 



4 












50 100 150 200 250 



A 



FIG. 8. The collisional width of the isoscalar giant quadrupole resonances obtained by using 
Eq. (40) with £(r, luo) from Fig. 7. The notations is the same as in Fig. 7. The experimental data 
were taken from Refs. [31, 32]. 
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FIG. 9. The diffusivity parameter a of the equilibrium distribution function in momentum 
space as a function of the distance r to the center of the nucleus for both harmonic oscillator 
(curve HO) and Woods-Saxon (curve WS) potentials. 
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FIG. 10. The equilibrium distribution function in the WS potential as obtained from the semi- 
classical expansion on the Hermite polynomials Eq. (59), with N = 3. 
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FIG. 11. The local damping parameter h/T2(r,u>o) for the nucleus with A = 224, calculated for 
a spherical WS potential. For curve 1 we use the equilibrium distribution function (59) and for 
curve 2 we use the Fermi distribution function (60). 
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FIG. 12. The damping factor £ in the case of the WS potential. The notations is the same as 
in Fig. 11. 
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FIG. 13. The collisional width of the isoscalar giant quadrupole resonances obtained by using 
Eq. (40) with £(r, ojq) from Fig. 12. The notations of curves 1 and 2 is the same as in Fig. 12. 
The dashed curve is from Fig. 8. The experimental data were taken from Refs. [31, 32]. 
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